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ABSTRACT 

We have observed the 5.54 s anomalous X-ray pulsar XTE J18 10-197 at radio, millimeter, and infrared (IR) 
wavelengths, with the aim of learning about its broad-band spectrum. At the IRAM 30 m telescope, we have 
detected the magnetar at = 88 and 144 GHz, the highest radio-frequency emission ever seen from a pulsar. At 
88 GHz we detected numerous individual pulses, with typical widths ^ 2 ms and peak flux densities up to 45 Jy. 
Together with nearly contemporaneous observations with the Parkes, Nanfay, and Green Bank telescopes, we 
find that in late 2006 July the spectral index of the pulsar was -0.5 < a ^ (with flux density S^, oc 
over the range 1.4-144 GHz. Nine dual-frequency Very Large Array and Australia Telescope Compact Array 
observations in 2006 May-September are consistent with this finding, while showing variability of a with time. 
We infer from the IRAM observations that XTE J1810-197 remains highly linearly polarized at millimeter 
wavelengths. Also, toward this pulsar, the transition frequency between strong and weak scattering in the 
interstellar medium may be near 50 GHz. At Gemini, we detected the pulsar at 2.2 /xm in 2006 September, at 
the faintest level yet observed, TiQ = 21.89±0.15. We have also analyzed four archival IR Very Large Telescope 
observations (two unpublished), finding that the brightness fluctuated within a factor of 2-3 over a span of 3 
years, unlike the monotonic decay of the X-ray flux. Thus, there is no correlation between IR and X-ray flux, 
and it remains uncertain whether there is any correlation between IR and radio flux. 

Subject headings: pulsars; individual (XTE J1810-197) — stars: neutron 



1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) are neutron stars with 
long spin periods and extremely strong surface magnetic 
fields, the decay of which is largely responsible for their 
observed X-ray emission according to the magnetar model 
(iDuncan & Tho mpson 1992). Of the dozen identified mag- 
netars (including soft-gamma repeaters in addition to AXPs), 
half have also been detecte d at infrared (IR) wavelengths (see 
IWoods & Thompsonll2006l for a review of magnetar proper- 
ties). While th e IR fluxes of mag netars have been observed to 
vary (e.g., HuUe man et alj|2004l) . it is neither clear how these 
variations relate to changes in observed X-ray flux, nor what 
is the origin in detail of the IR radiation. 

XTE J18 10-197 is a transie nt AXP with period 5.54 s, first 
detected (.Ibrahim et alj|2004b when its X-ray flux increased 
^ 100-fold co mpared to the historical le vel maintained for at 
least 24 years (iHalpern & Gotth elf 2005). Four years after the 
outbu rst, its X-ray flu x has decreased nearly to the quiescent 
level (iGotthelf & Hal pern 2007). It was also detected in the 
IR in 2003 October dlsrael et al. 2004), and was observed 5 
months later to be 60% fainter (Rea et al. 2004). 

XTE J 1 810-197 is the on ly magnetar known to emit radio 
waves (Halpern et alJl200"5h . This emission is entirely pulsed 
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(ICamilo et alJl2006l: iHelfand et alJl2007h. is al igned in phase 
with the X-ray pulsations ( Camilo et al.ll2007a ^. and is highly 



polarized, in some respects being similar to that of ordinary 
young radio pulsars (Camilo et al. 2007b; see also Kramer et 
al. 2007). However, unlike in ordinary pulsars, the radio emis- 
sion appears to be transient, the average flux density varies 
intrinsically by factors of up to ^ 3 on approximatel y daily 
times cales, and the pulse profiles vary as well (Camil oet alJ 
I2OO6'). By early 2007, the average radio flux density was 
a factor of ^ 20 lo wer than when pulsa tions were first ob- 
served 1 yr earlier (ICamilo et alJl2007ah . Also, unlike the 
vast majority of ordinary pulsars, XTE J1810-197 has a radio 
spectrum approximately consistent with being flat {a > -0.5; 
Si, oc i/"), o ver the large frequency range 0.7 < < 42 GHz 
(ICamilo et al. 2006). 

The radio window provides a new opportunity for learning 
about the emission mechanism(s) of magnetars. For instance, 
an extrapolation of the radio spectrum of XTE J181 0-197 ex- 
ceeds its IR flux, so that the radio and IR emission could have 
the same origin. However, the shortest radio wavelength at 
which a detection had been reported (7 mm) is 3000 times 
larger than the IR wavelengths, so that coverage at shorter 
radio wavelengths is desirable, also to provide further con- 
straints on the radio emission. Variability in the IR and com- 
parison with radio variability would also be of great interest, 
but only two IR observations have been reported, and they 
predate the discovery of radio pulsations. 

With these objectives in mind, we report here on observa- 
tions of XTE J18 10-197 at the IRAM telescope, which has 
the capability to detect pulsars in the 1-3 mm range, on si- 
multaneous radio observations, and also on a new IR observa- 
tion at Gemini, along with analysis of archival IR observations 
from the Very Large Telescope (VLT). 

2. OBSERVATIONS, ANALYSIS, AND RESULTS 
2.1. Millimeter wavelengths: IRAM 
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TABLE 1 

Sensitivity of XTE J 1 8 1 0-1 97 observations at IRAM 
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Fig. 1 .— XTE J1810-197 at IRAM at a central frequency of 88.5 GHz on 
2006 July 19 (left) and 144 GHz on July 20 (right). Only signal from one of 
the polarizations (H) is folded in these plots. The variation with time of the 
signal stre ngth (bottom plots) is likely mainly due to interstellar scintillation 
(see § 12. It . The white spaces correspond to times during which the telescope 
was not pointed at XTE J1810-197. and the pulse profiles displayed with 256 
phase bins are shown twice in each panel. 



We observed XTE J1810-197 at the IRAM 30 m telescope 
in Pico Veleta, Spain, during 2006 July 19-21''. We recorded 
four intermediate frequencies (IFs) simultaneously. These 
correspond to two linear polarizations (rotating with tele- 
scope elevation, but parallel and perpendicular to the ground 
at the inputs to the receivers, here designated H and V, respec- 
tively) at one of two possible pairs of frequencies (cryostats 
AB, with receivers operating at nominal frequencies of 100 
and 230 GHz, or CD, at 150 and 270 GHz). We used the 
AB receivers on July 19 and 21, and CD on July 20. We 
tuned all receivers in double side-band mode, with the IFs for 
the 100 GHz receivers centered at 1.5 GHz with 1 GHz band- 
width, and those for all others centered at 4 GHz with 2 GHz 
bandwidth. At the lowest frequency used, we set the local 
oscillator to 88.5 GHz, thereby recording combined signals 
from 87 ± 0.25 GHz and 90 ± 0.25 GHz. For the lowest fre- 
quency of the CD cryostat pair, signals were recorded from 
140 ±0.5 GHz and 148 ±0.5 GHz, for a central sky frequency 
of 144 GHz. The other central frequencies used were 224 and 
264 GHz. Dispersion is irrelevant at these frequencies, and 
we sampled the total power IF signals every 0.5 ms using a 
16-bit analog-to-digital converter, with the output recorded 
in VME modules using a custom configuration built by J. P. 
and colleagues at IRAM. The data stream containing the sam- 
pled power from the 4 IFs and a time stamp, referenced to 
the Observatory time standard and ultimately to UTC, was 

' IRAM operates both a six-element interferometer in France and the sin- 
gle dish in Spain; throughout we refer to the IRAM telescope as a shorthand 
for the latter. 



Date 
(MJD) 


(Jy) 


'^sys 

(Jy) 


'^sys 

(Jy) 


,,1^1 V 
'^sys 

(Jy) 


53935 


692-783 


775-843 


5640-7850 


6430-8960 


53936 


1780-2590 


1640-2420 


23000-53400 


23400-55000 


53937 


702-862 


775-942 


4290-10700 


4940-12100 


54079 


632-650 


725-729 


2130-2270 


2400-2560 


54118 


660-693 


757-815 


2650-2810 


2800-3000 


54172 


608-725 


706-832 


2250-2910 


2440-2800 


54180 


582-631 


687-740 


1900-2080 


2080-2290 


54207 


590-642 


698-750 


2700-3040 


2810-3120 



Note. — On MJD 53936, center frequencies were U] = 144 and 
U2 = 264 GHz, respectively; on all other days, 88 and 224 GHz. For each 
frequency, we provide the range of measured system equivalent flux den- 
sity s epara tely for the "horizontal" (H) and "vertical" (V) polarizations 
(see ji l2.lt . These were all corrected for air mass, measured hourly. The 
bandwidths used were 1 GHz at 88 GHz and 2 GHz at other frequencies. 




53936 53937 
Date (MJD) 

Fig. 2. — Peiiod-averaged flux densities of XTE J1810-197 across the fre- 
que ncy range 1.4-144 GHz obtained over a period of 2.5 days (see §§ 12. H and 
12.2) , plotted as a function of frequency (top) and date (bottom). For each com- 
bination of frequency and telescope we use a unique symbol that is common 
to both plots. At 1.4 GHz, multiple symbols represent the following data: 
Parkes (open triangles); Nan^ay (solid triangle); VLA (open square). With 
the exception of the VLA measurements whose flux density uncertainties are 
recorded in Table|2] the fractional uncertainties are ~ 25%. 



recorded to disk for off-line analysis. A typical observing 
session of 5 hr consisted of sets of pointing, focus, and flux- 
calibration scans (on Mars, Jupiter, or Venus), interspersed 
with four scans on XTE J18 10-197, each about 1 hr in length, 
at elevations of 12° < e < 39°. Weather conditions were vari- 
able, and the sensitivity varied on each night by up to 20% at 
88 GHz and 50% at 144 GHz. 

Due to large changes in the atmospheric conditions at 
IRAM on timescales of tens of seconds to several minutes, 
folding each XTE J1810-197 time-series modulo the 5.54s 
pulse period (using the contemporary radio ephemeris; § |2.2t 
resulted in a strongly varying off -pulse baseline which greatly 
reduced the sensitivity to the detection of pulsations. In ad- 
dition, periodic interference contaminated the H polarization 
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TABLE 2 

Radio and millimeter observations of XTE J1810-197 



Date 
(MJD/yymmdd) 


Time 
(hr) 


Frequency v 
(GHz) 


(mJy) 


Telescope 


53862.4/060507 


0.2 


1.4 


9.8 ±0.4 


VLA 


53862.4/060507 


0.2 


8.5 


5.7 ±0.2 


VLA 


53873.3/060518 


0.6 


1.4 


3.3±0.2 


VLA 


53873.3/060518 


0.6 


4.9 


3.5 ±0.3 


VLA 


53889.4/060603 


0.7 


1.4 


8.1 ±0.5 


VLA 


53889.4/060603 


0.7 


4.9 


6.2 ±0.2 


VLA 


53894.6/060608 


7 


1.4 


8.5 ±0.2 


ATCA 


53894.6/060608 


7 


2.4 


8.5 ±0.2 


ATCA 


53897.5/060611 


0.7 


1.4 


10.4 ±0.5 


VLA 


53897.5/060611 


0.7 


4.9 


7.5 ±0.2 


VLA 


53905.5/060619 


0.3 


1.4 


8.0 ±0.9 


VLA 


53905.5/060619 


0.3 


4.9 


8.0 ±0.2 


VLA 


53906.1/060620 


0.9 


0.35 


9.4 


GBT 


53935.5/060719 


0.3 


1.4 


4.9 


Parkes 


53935.5/060719 


0.6 


1.4 


2.6 


Parkes 


53935.6/060719 


2.3 


6.6 


0.7 


Parkes 


53935.7/060719 


0.4 


1.4 


1.9 


Parkes 


53935.9/060719 


0.2 


1.4 


1.8 


Nanfay 


53935.9/060719 


3.6 


88.5 


1.2 


IRAM 


53936.0/060720 


1.5 


9.0 


0.7 


GBT 


53936.1/060720 


0.9 


19.0 


1.0 


GBT 


53936.3/060720 


0.7 


1.4 


2.0 ±1.0 


VLA 


53936.3/060720 


0.7 


4.9 


1.4±0.2 


VLA 


53936.9/060720 


4.7 


144.0 


1.2 


IRAM 


53937.9/060721 


1.7 


88.5 


1.2 


IRAM 


53958.2/060811 


0.7 


1.4 


6.0 ±0.6 


VLA 


53958.2/060811 


0.7 


4.9 


2.7 ±0.2 


VLA 


53983.2/060905 


0.7 


1.4 


3.4 ±0.6 


VLA 


53983.2/060905 


0.7 


4.9 


1.1 ±0.3 


VLA 


53991.8/060913 


0.3 


1.4 


1.6 


Nangay 


53992.7/060914 


0.2 


1.4 


1.6 


Nanfay 


54069.5/061130 


1.1 


1.4 


0.4 


Nanfay 


54069.9/061130 


0.6 


0.35 


1.6 


GBT 


54078.5/061209 


0.7 


1.4 


0.7 


Nanfay 


54079.5/061210 


3.6 


88.5 


0.22 ±0.11 


IRAM 


54118.4/070118 


4.0 


88.5 


<0.1 


IRAM 


54172.2/070313 


4.5 


88.5 


<0.1 


IRAM 


54180.2/070321 


4.0 


88.5 


<0.1 


IRAM 


54207.2/070417 


3.2 


88.5 


0.34±0.17 


IRAM 


54208.2/070418 


0.8 


1.4 


1.2 


Nanjay 



Note. — We list here all our simultaneous dual-frequency obser- 
vations, along with selected others, ordered by date. At IRAM, the fre- 
quency pairs were 88.5 and 224 GHz, or 144 and 264 GHz . Th e pulsar 
was detected at only the lower frequency of each pair (see 8 12.11 . Unless 
specifically given, the fractional flux density uncertainties are about 25%. 





0.04 0.06 
Pulse Phase 

Fig. 3.— Single pulses of XTE J1810-197 at 88 GHz detected at IRAM on 
2006 July 19. The 10 single pulses, displayed with 0.5 ms resolution, were 
chosen from among the brightest in our data and are time-ordered (though not 
consecutive) from top to bottom. At bottom is the average pulse profile on 
this day, displayed with 1 1 ms resolution (see also Fig.[T). Only 10% of pulse 
phase is shown in order to emphasize the shaipness of the sub-pulse structure 
(with typical width ~ 2 ms); no emission was detected from other phases. 
Emission such as seen in the top five pulses trailing the main component(s) 
by ~ 0.07 in phase was among the brightest in the dataset, but was infrequent 
enough that it did not contribute significantly to the average profile. The 
brightest pulse detected is shown fourth from the top, with a signal-to-noise 
ratio of 45 coiresponding to a peak flux density of approximately 45 Jy. As 
in Figs. [T] and |4] only signal from one polarization is displayed here. The 
high-frequency noise readily apparent in the baseline of the profiles is due to 
imperfectly excised interference. 



and severely degraded the V polarization, greatly increasing 
the off-pulse noise levels in the folded profiles. We believe 
that these signals, at 50 and 100 Hz, as well as many harmon- 
ics of 1 Hz, are locally-generated compressor-related interfer- 
ence. 

In order to combat these issues, we filtered the data in two 
ways. First, to remove the low-frequency noise below ^ 
0.1 Hz (due primarily to the changing atmosphere), we high- 
pass-filtered each time series using a 3rd-order Bessel filter. 
We then Fourier transformed the time series and clipped the 
Fourier amplitudes of the strong 1 Hz harmonics as well as the 
prominent 50 and 100 Hz signals. Finally, we inverse Fourier 
transformed the data to regenerate mostly interference-free 
time series, which we then folded. After filtering, the V po- 
larization sensitivity to pulsed signals remained significantly 
worse than that of the H polarization. The pulsar was not de- 
tected in the V polarization. In the H polarization, the pulsar 
was clearly visible on both days at 88 GHz (3.4mm) and at 
144 GHz (2.1 mm; Fig.[T]), and was not detected at higher fre- 
quencies. 

We determined the period-averaged flux density by mea- 



suring the area under the profiles and dividing by the pulsar 
period. This was converted to an absolute Jansky scale us- 
ing the hourly flux-calibration scans, from which the system 
equivalent flux density 5sys was calculated for each frequency 
and polarization. The off -pulse profile rms then corresponds 
to Ssys(BT)~^/^ for a bandwidth B and an integration time T, 
where Ssys is corrected for air-mass attenuation exp(-T/ sin e), 
and T is the measured zenith opacity. The ranges of ^sys ob- 
tained for each observation are listed in Table [T] 

It is clear from Figure [T] that the observed flux density can 
vary greatly with time. Changing weather conditions and el- 
evation account for only up to 10% of this variation in the 
July 19 data, displayed in the left panel of the figure (see Ta- 
ble [T]). A great portion of the variation is most likely caused 
by interstellar scintillation (see iCamilo et al.u2006) . As an 
average, we estimate the flux density to have been 1.2 mJy 
at both 88 and 144 GHz, with a fractional uncertainty on the 
absolute values of about 25% (see Fig. |2]and Table [2]i. In 
order to calculate approximate flux density limits for non- 
detections, we assume a threshold signal-to-noise ratio of 5 
and a pulse duty cycle of 2%. For example, on July 19 we 
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Parkes, 53935.5, 1.4 GHz, (288 MHz, 0.27 h) 



Parkes, 53935.6, 6.6 GHz, (576 MHz, 2.3 h) 
GBT, 53936.0, 9 GHz, (800 MHz, 1.4 h) 

A 

GBT, 53936,1, 19 GHz, (800 MHz, 0.95 h) 
IRAM, 53935.8, 88 GHz, (1000 MHz, 1.1 h) 
RAM, 53936.8, 144 GHz, (2000 MHz, 3.6 h) 



Pulse Phase 

Fig. 4. — Average pulse profiles of XTE J1810-197 obtained at frequen- 
cies spanning 1.4— 144 GHz over a period of 1.5 days, at Parkes, GBT, and 
IRAM. Each profile is labeled by the telescope used, the date (MJD) of the 
observation, the central frequency and, in parentheses, the bandwidth and in- 
tegrat ion t ime used. IRAM profiles include data from only one polarization 
(see § 12.11 , and are displayed with 256 phase bins. Parkes and GBT profiles 
have 1024 bins. Profiles were aligned by eye such that the peak pulse com- 
ponent arrives near phase 0.15. The FWHM of these profiles are, in order of 
increasing frequency, 4.2, 2.2, 1.7, 1.6, 1.3, and 1.3% of the pulse period. 



obtain ^224 ^0.9 mJy. 

The observed flux density was so great on July 19 toward 
the end of the third scan (Fig. [T]) that we were able to detect 
numerous individual pulses from XTE J1810-197 at 88 GHz. 
On this day we detected single pulses with signal-to-noise ra- 
tio > 4 from about 15% of all rotations of the neutron star 
(Fig. [3]). The largest pulses had a peak flux density of ^ 45 Jy, 
comparable to the strongest celestial sources known at 3 mm 
with the exception of the Sun, Jupiter, and Venus (although 
only for ^ 1 ms out of every ^ 5 s). On July 21 we also de- 
tected single pulses at 88 GHz. 

Following these detections, we monitored the pulsar at 
88/224 GHz on five occasions between 2006 December and 
2007 April, using identical observing parameters. These lat- 
ter observations were done in winter and early spring, under 
better weather conditions (Table [U. The pulsar was detected 
on two of these occasions, at 88 GHz in the H polarization, in 
December and April, with S^s = 0.2-0.3 mJy (see Table|2|. In 
2007 April, 5224 < 0.4 mJy. 

2.2. Radio: Parkes, Nangay, GBT, VLA, andATCA 

Because the intrinsic flux density and pulse profile of 
XTE J1810-197 vary on 1 day timescales, determina- 
tion of a spectral index ideally requires simultaneous multi- 
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Date (MJD) 



Fig. 5. — Spectral indices for XTE J1810-197 determined from simul- 
taneous dual-frequency VLA (triangles) or ATCA (scjuare) observations, or 
near-simultaneous GBT data (crosses). For each simultaneous observation 
the lowest frequency used was 1 .4 GHz. For the VLA the second frequency 
was 8.5 GHz for the first point plotted, and 4.9 GHz for all others (see Ta- 
ble |2). The point with the largest error bar was obtained on 2006 July 20. 
The first GBT point is based on flux density data obtained at 1.9 and 9 GHz 
over a 4hr interval. The last point is based on observations at 0.8, 1.9, 5, and 
9 GHz, obtained over 7 hr Error bars are 1 a. 



frequency observations. We therefore observed at Parkes, 
Nan9ay, the Green Bank Telescope (GBT), and the Very Large 
Array (VLA) on 2006 July 19-20, nearly simultaneously with 
IRAM. Table |2] lists these and all dual-frequency observa- 
tions at the VLA and the Australia Telescope Compact Array 
(ATCA). The methods used have been described elsewhere 
(see Camilo et al. 2006 for Parkes and GBT; Camilo et al. 
2007a for Nan9ay and VLA; Helfand et al. 2007 for ATCA). 

The 1 .4 GHz flux density decreased over a period of a few 
hours shortly before the first IRAM observation (Fig. |2] and 
Table|2]l. The average pulse profiles during these different flux 
states were identical (we show the first one in Fig. |4|i, unlike 
an observation in 2006 S eptember w ith sudden simultaneous 
flux and profile changes jCamilo et al. 2007a). 

The nominal flux density of the pulsar at 19 GHz was 
slightly greater than at 9 GHz (see Fig. |2]l. However, the ob- 
served flux at ly > 9 GHz varies on short timescales due to 
interstellar scintillation (Camilo et al. 2006; Fig. [T]), and this 
may bias somewhat those flux densities obtained here (i.e., 
1 hr does not correspond to many scintles). 

Using simultaneous 1.4 and 4.9 GHz observations on 2006 
July 20, we measured a spectral index a = -0.3 ± 0.4. The 
large uncertainty is due to the uncertain 1.4 GHz flux den- 
sity, the smallest ever measured at the VLA for the magnetar. 
On a total of nine occasions we have used accurately flux- 
calibrated simultaneous dual-frequency data to determine the 
spectral index of XTE J1810-197, shown in Figure |5] From 
this, it is apparent that the radio spectrum, averaged over all 
its pulse profile components, varies, and there is a hint that 
since late 2006 July it might have become steeper (at least in 
the 1.4-4. 9 GHz range). This appears to be supported by esti- 
mates of a obtained from near-simultaneous multi-frequency 
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TABLE 3 



Infrared 


OBSERVATIONS 


OF XTE J1810- 


197 


Date 


Exposure 


magnitude 


Kels. 


(MJD/yymmdd) 


(minutes) 






52920.05/031008 


32 


20.9 ±0.2 


(1) 






20.8 ±0.1 


(2) 


53078.37/040314 


36 


21.21 ±0.14 


(1) 






21.36±0.07 


(3) 


53258.99/040910 


36 


20.9 ±0.2 


(1) 


53259.99/040911 


60 


20.82 ±0.16 


(1) 


53992.21/060914 


38 


21.89±0.15 


(1) 



Refer ences. — (1) This work; (2) [Israel et ai] f2004l) ; (3) 
IRea etalj 12004) . 



Note. — Two entries are provided for each of the first two ob- 
servations, con'esponding respectively to our analysis, and to the 
original published analysis. The first four observations listed were 
obtained at the VLT, while Gemini was used for the last one. Dates 
correspond to observation start times. All observations used adap- 
tive o ptics , and the magnitudes listed are on the 2MASS system 
(see § 12.3) . Uncertainties are all given at the 1 a level. 

pulsed flux measurements at GBT: those shown in Figure|5]for 
two epochs in 2006 October and November hav e a < -0.5, 
while before mid-2006, a > -0.5 (ICamilo et al. r2006). We 
have also extended the range of frequencies over which the 
magnetar has been detected, with two GBT observations at 
0.35 GHz. Together with observations at other frequencies 
within half a day of these (Table |2|i, we obtain a = in 2006 
June, but a = -l in 2006 November. This further supports the 
notion that the spectrum of XTE J1810-197, while variable 
with a ranging between approximately and -1, may have 
become generally steeper after mid-2006. 

2.3. Infrared: Gemini and VLT 

We obtained a near-IR ^Tj-band (2.15 /im) observation of 
XTE J1810-197 at Gemini-North on 2006 September 14. 
We used the adaptive optics (AO) system Altair with the 
near-IR imager NIRI (lHodappetal.ll2003h . With this con- 
figuration, the 1024 X 1024 pixel Aladdin InSb array cov- 
ers 22 X 22arcsec^ at 21.9maspixer'. For photometric cal- 
ibration, as well as the astrometric analysis described in 
iHelfand eT al. (2007), we also analyzed /Tj-band observations 
taken on 2003 September 18 with NIRI on Gemini without 
the AO system (for which the detector covers 2x2 arcmin^ at 
1 17 mas pixel"'; 19minute exposure time). 

For comparison, we analyzed observations taken at the 
VLT using the AO system NAOS with the CONICA camera 
dlenzen et all 120031: iRousset et al.1 l2003h . NAOS-CONICA 
also uses a 1024 x 1024 Aladdin detector, covering 27 x 
27arcsec^ at 27.0 mas pixe l"'. These o bservations were ob- 
tained in 2003 October (llsrael et all [2004 ). 2004 March 
(iRea et al.i r2004). and 2004 September (on two nights, pre- 
viously unpublished). We summarize in Table[3]the IR obser- 
vations of XTE J1810-197 analyzed here in detail. 

All observations were taken in a similar way, with images 
taken at dithered positions, and each image consisting of one 
or more co-added exposures, with the counts in each expo- 
sure determined from the difference between series of read- 
outs before and after the actual integrations. We reduced all 
observations in an identical fashion, using the Munich Image 
Data Analysis System (MIDAS) or the Image Reduction and 
Analysis Facility (IRAF). We corrected for pixel-to-pixel sen- 
sitivity variations using sky flats constructed from the images 
themselves, aligned the images to integer pixel boundaries, 
and took averages. The 2003 October data had unrepealed 
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Fig. 6. — Adaptive optics images of XTE Jl 8 10-197 at 2.2 ^m taken with 
the VLT in 2004 March (left) and with Gemini in 2006 September (right). 
The tickmarks are centered on the accurate VLB A po sition. Between the two 
epochs, the infrared counterpart has faded (see j) |2.3| for details). 



horizontal "banding" that was not removed by the flat-fielding 
process. We removed this feature using the median value 
for ea ch data row before stacking images. In IHelfand etal] 
("2007), we used the 2004 and 2006 observations in Table[3]to 
confirm that the IR counterpart has the same proper motion as 
that measured from radio VLBA observations. 

In Figure |6] we compare the final VLT image from 2004 
March with our new Gemini image from 2006 September 
Over this 2.5 yr period, the IR counterpart has moved about 
one pixel to the south relative to local stars due to its proper 
motion, and it has faded. For a quantitative an alysis of all th e 
images, we used the DAOPHOT II package (IStetsonlll987l) . 
running inside MIDAS and IRAF, to derive instrumental mag- 
nitudes by fitting a model point-spread function (PSF) for 
stars on the average im ages. Following the recommenda- 
tions of IStetsoiil (Il987b . the PSF was derived from brighter 
stars in an iterative fashion, where in each iteration nearby 
fainter stars are removed with the current best model. For a 
well-matched model PSF, with minimal residuals, we chose 
a Lorentzian analytical base, and linear dependence of the 
shape on position. We determined uncertainties on the instru- 
mental magnitudes by adding fake stars of a range of bright- 
ness to the images, re-measuring those, and determining the 
rms deviations. We found that these numbers were very simi- 
lar to the measurement uncertainties returned by DAOPHOT. 

We calibrated our instrumental magnitudes relative to 
the 2MASS (Two Micron All Sky Survey) catalog 
(ISkrutskie et al .1120061) in two steps. First, we derived offsets 
for 2MASS stars measured in the wider 2003 September NIRI 
image. We found that a constant offset gave a good fit, with 
the rms residual of 0.06 mag being only slightly larger than the 
typical 2MASS measurement error of typically ~ 0.05 mag 
for the 21 stars we used. Thus, the photometry of this image 
is tied to the 2MASS system to about 0.01 mag. 

We then transferred the photometry to the AO images using 
fainter stars. For this step, we found that for most images the 
residuals around the mean of instrumental minus calibrated 
magnitudes had rather large rms, of 0.02, 0.08, 0.06, 0.06, 
and 0.10 mag respectively for each the five data sets (see Ta- 
ble O, with no obvious dependencies on brightness, position, 
or proximity to the guide star used. Since most of these rms 
residuals are much larger than the measurement uncertainties, 
we use them as a measure of the calibration uncertainty for 
any given star Likely, the large scatter reflects the general 
difficulty of doing reliable photometry on AO-corrected im- 
ages. Indeed, the largest scatter is found for the 2004 March 
VLT and the 2006 Gemini data, which had the best AO cor- 
rection, while the smallest is for the 2003 October VLT data. 



6 



CAMILO ET AL. 



which had the worst correction. 

The Ks magnitudes thus obtained for XTE J1810-197 are 
Hsted in Table |3] with measurement and caHbration uncer- 
tainties added in quadrature (all uncertainties here are given 
at the 1 a confidence level). The values for 2003 Octo ber and 
2004 M arc h are consistent w ith those reported by Israel et alJ 
(120041) and iRea et all ( |2004|) . respectively, but our uncertain- 
ties are larger by a factor of 2. This cannot be due solely 
to the fact that we included the calibration uncertainty, since 
for the two observations concerned this is smaller than the 
measurement uncertainty for the faint target. In order to ver- 
ify our estimates of the uncertainties, we checked the mea- 
surements of three other, relatively faint stars (at « 2.0" SSE, 
1.7" ENE, and 2.8" N of the target; the former two are vis- 
ible in Fig. |6l). For the five observations, their magnitudes 
are [20.35(12), 20.39(11), 20.30(12), 20.33(11), 20.39(10)], 
[20.30(11), 20.33(11), 20.23(11), 20.27(10), 20.40(10)], and 
[21.5(4), 20.94(14), 21.0(2), 21.0(2), 21.01(11)], where the 
uncertainties again include the calibration uncertainties. The 
uncertainty estimates are comparable to the scatter between 
the magnitudes, suggesting that they are reasonable (if per- 
haps slightly high). 

Independent of the details of the uncertainties, our measure- 
ments show that the IR brightness of XTE J1810-197 did not 
vary greatly in 2003 and 2004, by < 50%, declining more 
sharply only by 2006. Overall, these (sparsely sampled) ob- 
servations do not show a clear trend. We also searched for 
intra-night variations in the 2004 March data, by considering 
the observation in two separate 18 minute segments, finding a 
magnitude difference of 0. 1 1 ± 0.23. At the ^ 1 hr or ~ 1 day 
timescales, we thus have found no evidence for variations at 
the > 20% level. 

Fina lly, using the /Tj-band zero point of ICohen et alJ 
(l2003h . the Gemini measurement in 2006 September corre- 
sponds to a flux density of 1.17 ± 0.17 /iJy. If A y = 3.6 
(iHalpern & Gotthelf 2005), then Ak = 0.40 (Schleg el et aP 
Il998l) . so that the de-reddened flux is 1.45 times larger, 
1.69±0.25AiJy. 

3. DISCUSSION 

The remarkable detection of XTE Jl 8 10-197 at frequen- 
cies of 88 and 144 GHz wit h IRAM is a record among pul- 
sars (see iMorris etal.lll997L for previous highest-frequency 
detections). This magnetar has now been detected at radio- 
millimeter wavelengths spanning a factor of 400. The spec- 
trum over this expanse is flat compared to that of most pul- 
sars (a_w_^L6jmeasured over a narrower frequency range; 
see iLorimer et alj|1995h . but it is difficult to make a defini- 
tive quantitative sta tement: (1) as for ordinary pulsars (e.g., 
iKramer et alj2003l) . individual phase components of the pulse 
profile may have differing spectra; (2) at high frequencies 
(> 9 GHz) there is considerable variability in measured flux 
density due to interstellar scintillation; (3) the intrinsic flux 
varies with time; (4) the spectrum also apparently varies with 
time (see Fig. |5]i; and (5) the intrinsic spectrum need not be 
represented by a single power law. 

The only simultaneous multi-frequency detections among 
those summarized in Figure |2] were the ones at the VLA (1.4 
and 4.9 GHz), and 1.4 GHz at Nan9ay along with 88 GHz 
at IRAM. The spectral index obtained from the first pair is 
a = -0.3, while from the second, a = -0.1, in both cases with 
substantial uncertainties and also suffering from some of the 
problems noted above. Attempts to infer a from other mea- 
surements collected within a period of 2.5 days (see Fig. |2]i. 



also suffered from such problems. (See also Fig. |5]) With 
the foregoing caveats in mind, we summarize the situation 
thus: before 2006 August, the pulse-averaged spectrum of 
XTE J1810-197 over the range 1.4-144GHz could be use- 
fully described by a single spectral index in the range -0.5 < 
a < 0, with some time variability observed. Since then, the 
spectrum has apparently steepened, with -1 -0 < a < -0.5. 

In retrospect, the flux of the pulsar started changing dramat- 
ically around mid 2006 July, near the time of the first IRAM 
observations. This was accompanied by a change in the char- 
acter of pu lse profile variations , as well as large torque vari- 
ations (see lCamilo et al.|[2007al) . The flux density at 1.4 GHz 
by early 2007 was generally quite low (< 0.5 mJy — although 
fluctuations continue: see the last entry in Table|2]l. Together 
with a moderate spectral index of a < -0.3, this can account 
for the non-detections at > 88 GHz in 2007 (Table |2l). 

The average pulse profiles of XTE J1810-197 observed 
during 2006 July 19-20 and spanning 1.4-144 GHz are shown 
in Figure|4] The profiles at 6.6-19 GHz are very similar, while 
in all 1 .4 GHz observations on this date the "precursor" com- 
ponent is absent (phase 0.85 in the figure). This is unexpected, 
because the Parkes 1.4 GHz observations bracketed in time 
that at 6.6 GHz, ensuring that the differences observed be- 
tween the pulse profiles were not caused by a global change 
in the profile. Rather, the absence of the precursor at 1 .4 GHz 
may point on this day either to an extremely positive spec- 
tral index for that component (a > 1), or to absorption of the 
1.4 GHz radiation along the line of sight at the location in the 
magnetosphere corresponding to this pulse phase. 

The absence of the precursor at 88 GHz may be explained 
by a spectrum over 19-88 GHz steeper than that of the main 
component by Aa < -0.7. In making this estimate we as- 
sumed that roughly half of the precursor power would be 
present in each of the H and V polarizations. This follows 
from a computation of the polarimetric Stokes parameters as 
they would appear for the IRAM observations, under the as- 
sumption that the precursor was polariz ed at millimeter wave - 
lengths as it was earlier at 1.4 GHz (see lCamilo et al.|[2007bl) . 
A similar computation for the main profile component using 
the (different) absolute position angle of linear polarization 
observed at 1 .4-8.4 GHz (Camilo et al. 2007b), shows that the 
vast majority of its power in the IRAM observations should 
be present in the H polarization, as is indeed the case (§ 12. 11 ). 
This provides evidence that XTE J1810-197 remains highly 
linearly polarized at wavelengths as short as 2 mm. For com- 
parison, some ordinary pulsars appear to depolarize signifi- 
cantly with decreasin g wavelength, down to the limit of the 
observations at 1cm jXilouris et al.l [19961) . However, sep- 
arate p rofile components can ev olve differently with wave- 
length. iKarastergiou et al.] (12005 ^ consider a model where or- 
thogonal polarized modes have different spectra, leading nat- 
urally to some highly polarized components with flat spectra, 
down to 10 cm wavelength. Also, several young pulsars re- 
main highl y polarized down to the limit of the observations 
at 3.5 cm (J ohnston et alj|2006l) . It is therefore possible that 
in this respect XTE J1810-197 differs from ordinary young 
pulsars more in having a flat spectrum all the way down to 
2 mm, which enables us to make the measurements, than in 
remaining highly p olarized at those wavelengths (see also 
ICamilo et al.ll2007bl) . In any case, these two features may not 
be independent. 

As seen clearly in Figure[Tl the flux density of XTE J1810- 
197 observed at millimeter wavelengths through a bandwidth 
of 1-2 GHz varies by factors of a few on timescales of 1 hr 
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Most of this variation is likely caused by interstellar scin- 
tillation (§ 12.1b . The surprise is not that the received flux 
varies, but that it varies so relatively little. Qualitatively, the 
flux modulation for these IRAM observations is much smaller 
than that observed at 42 G Hz and more com parable to that 
observed at 14-19 GHz (cf. lCamilo eTaHIIOOe) . Perhaps for 
this object, the transition between strong and weak interstel- 
lar scatterin g (at which the flux modulation peaks; see, e.g., 
iLorimer &" Kramer 2005) is about 50 GHz. Detailed analysis 
of scintillating behavior for XTE J18 10-197 will be presented 
by Ransom et al. (in preparation). 

The single pulses seen from XTE J1810-197 at 88 GHz on 
2006 July 19 (Fig.|3]l appear qualitatively simila r to those ob- 
serve d at 2 and 42 GHz about 2 months earlier dCamilo et alj 
12006). In particular, they are narrow (a few ms), wander about 
with no obvious phase coherence from rotation to rotation, 
and gradually build up the much wider average profile. Also, 
some particularly strong pulses are emitted from rotational 
phases at which the average pulse is weak or not detectable, 
implying significantly different pulse energy distributions at 
different phases. A detailed treatment of single-pulse behav- 
ior from XTE J1810-197 will be presented elsewhere. 

With the detection of XTE J18 10-197 down to 2 mm wave- 
length, and shortward of 2.2 /im, there remains a factor of 
1000 in wavelength where it is undetected (Wang et al. 2007 
report upper limits at 24, 8, and 4.5 /im). As discussed earlier, 
the radio-to-millimeter spectrum is somewhat uncertain, and 
is variable, but it is flat enough that the radio spectrum extrap- 
olated to IR wavelengths exceeds the IR fluxes (possibly by a 
very large amount; see Fig. |7]l. Therefore, it is possible that 
the radio spectrum steepens smoothly to join the IR. Contem- 
poraneous detections in several IR bands would be useful to 
delineate the shape of the spectrum in this region. Similarly, 
detection of IR pulsations from XTE J1810-197 would be of 
great help in pinning down the emission mechanism, but this 
has not been attempted due to the faintness of the source. 

The X-ray spectra in Figure|7]represent blackbody fits to the 
high and low states observed histor ically, discussed in more 
detail bv lHalpern & Gotthell) (l2005l) . These are corrected for 
interstellar absorption, which has been fitted as a free param- 
eter to the X-ray spectrum. The high-state spectrum requires 
two blackbodies to get a satisfactory fit, while detailed mod- 
els involving the effects of the stellar atmosphere and reso- 
nant cyclotron scattering in the magnetosphere can reproduce 
the spectrum using one surface temperature and one magne- 
tospheric electron temperature (Giiveret al. 2006, 2007). In 
either case, it is evident that the Rayleigh-Jeans tail of the 
thermal X-ray spectrum cannot account for the IR flux, which 
must come from a separate mechanism. (Until recently, it was 
customary to fit X-ray spectra of AXPs with a steep power- 
law plus a blackbody, but this is unphysical because it grossly 
overpredicts the optical and IR fluxes, and because Comp- 
tonization of X-ray blackbody photons can only produce high- 
e nergy powe r-law tails, not low-energy ones.) 

iRea et alj (12004) noted that the IR and X-ray flux from 
XTE J1810-197 had both decreased by a factor of about 2 
within a period of 5 months (see the first two observations in 
Table|3j. They considered that a fossil disk reprocessing some 
X-rays might account for this correlation. However, we have 
shown that the IR emission from XTE J1810-197 fluctuated, 
rather than following a monotonic decreasing trend (§ 12. 3t . In 
2006 Septembe r, it wa s a factor of about 2 fainter than when 
observed bv Rea et al.l ([2004.) 2.5 years before (Table |3]l. For 
comparison, the X-ray flux decreased monotonically by a fac- 
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spin-down luminosity of XTE J18I0-197 in early 2007 was approximately 
2 X 10" ergss-' (see lCamilo et aU2007iO . 



tor of ~ 20 in the same period (iGotflielf & HalpernI 12007'). 
Therefore, the IR and X-ray fluxes are not simply correlated. 

Four AXPs, including XTE J1810-197, have displayed IR 
variations, although no general trend is evident thus far In 
IE 2259H-586, following an X-ray outburst and accompa- 
nying IR increase, both decayed simil arly and the IR flux 
reached the quiescent level after 1 yr (iTam et al.ll2004l) . In 
IE 1048.1-5937, large IR variations are not correlated with 
the spin-down rate, but may be anticorrelated with the X-ray 
flux ©urant & van Kerkwijk 2005). And in 4U 0142H-61, al- 
though the IR, optical, and X-ray fluxes (and spectra) all vary, 
including IR variations of over a mag nitude on a timescale of 
days, there are no clear correlations (iDurant & van KerkwiikI 
2006). In XTE J 1 8 1 0-1 97, the latest IR flux is also the faintest 
yet observed, and the p resently observed sp in-down rate is 
the smallest on record dCamilo et alj|2007al) . The previous 
IR observations, brighter and perhaps showing some variabil- 
ity (Table Is), all took place within 1 .8 yr of the original large 
X-ray outburst ( Ibrahim et al.ll2004b . It is therefore possible 
to postulate that the latest IR detection corresponds to some 
sort of quiescent le vel, much like the st ar has approximately 
reached in X-rays dGotthelf & Halpernll2007) . while the ear- 
lier observations reflected at least in part transient responses 
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from an active magnetar, related to the outburst that led in 
the first place to the discovery of XTE J1810-197. However, 
the sampling of the IR observations is limited, and no clear 
conclusion can yet be reached. 

Likewise, we do not yet know whether the IR and ra- 
dio fluxes are correlated. The original radio detection, in 
2004 January, had a flux density at 1.4 GHz of 4.5 mJy 
dHalpern et al.ll20()5l) : at the time of the latest IR observation 
in 2006 September, this was near 1.6 mJy (Table |2]i; and be- 
tween 2006 February and September, the flux density ranged 
between 1 3 and 1 mJy, with day-to-day fluctu ations by fac- 
tors of up to 3 and a general downward trend (ICamilo et al.l 
l2007ah . At least in 2006, therefore, the radio flux fluctuated 
by a greater factor than did the IR between 2003 and 2006. 
But given the limited IR and nearly non-existent simultaneous 
radio-IR sampling, it remains possible that the fluxes in these 
two bands may be related. Simultaneous IR and radio ob- 
servations sampling a variety of timescales should settle this 
question. 
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